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ABSTRACT 


Two  new  transient  tasting  tschniquts  were  evaluated;  the  cantroid 
method  davalopad  by  Kohlmayr  and  the  cine  aero  interest l  technique.  The 
taro  intercept  aethod  was  found  to  be  the  aost  promising  of  the  two  but 
is  limited  to  values  of  Ntu  <  2.5.  The  centroid  t  ,-chnique  can  be  used 
effectively  when  the  value  of  H  is  less  than  5.0. 

A  heater  systoa  aade  of  .001  inch  dlaaeter  nlchrome  wire  was  designed 
and  tested  to  determine  its  effect  on  the  transient  testing  of  matrix 
type  heat  exchangers.  Because  the  design  showed  no  improvement  in  the 
test  results  and  was  unreliable  its  use  was  discontinued. 
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English  Letter  Symbols 


NOMENCLATURE 


* 

S 

b 


C 

8 


C 

P 


c 

8 


E 

G 


g 


h 


1(g) 

k 

k 


Matrix  total  heat  transfer  area 

Matrix  uni  ’<  -i  free  flow  area 

Matrix  total  fiontal  area 

Matrix  solid  cross-sectional  area 
available  for  thermal  conduction 

Fin  thickness 

Flow  passage  perimeter  (  fS  A^) 

Fluid  stream  thermal  capacity  rate 
(mcf) 

Fluid  specific  heat 

Matrix  thermal  capacity  (W^c^) 

Fluid  specific  heat  at  constant 
pressure 

Matrix  material  specific  heat 
Flow  passage  hydraulic  diameter 

(*rh) 

Friction  power  per  unit  area 

Flow  stream  mass  velocity  (m/A  ) 

c 

Proportionality  factor  in  Newton's 
Second  Law 

Normalized  fluid  temperature  at 
inlet  to  test  section 

Surface  heat  transfer  coefficient 
for  convection;  heat  transfer 
power  per  unit  area  per  degree 
temperature  difference 

Deviation  from  step 

Fluid  thermal  conductivity 

Matrix  thermal  conductivity 


sq  ft 
sq  ft 
sq  ft 
sq  ft 

ft 

ft 

Btu/(hr  deg  F) 

Btu/(lbo  deg  F) 

Btu/deg  F 
Btu/ (lbm  deg  F) 

Btu/(lbra  deg  F) 
ft 

hp/sq  ft 

lbm/ (hr  sq  ft) 

32.2  (lbm  ft)/(lbf 
sec*) 

dimensionless 
Btu/ (hr  sq  ft  deg  F) 

dimensionless 

Btu/ (hi  sq  ft  deg  F/ft) 

Btu/(hr  sq  ft  deg  F/ft) 


9 


L  Total  matrix  flow  length 

m  Has*  flow  rate 

P  Pressure 

p  Matrix  porosity  (A  /A,  } 

c  i  r 

q  Heat  transfer  rate 

R  Gas  constant  (53.35-alr) 

r^  Hydraulic  radius  (A£L/A) 

t  Temperature 

t*  Normalised  temperature 

u  Flow  velocity 

V  Matrix  volume 

m 

W  Matrix  mass 

s 

Fluid  mass  in  matrix 

x  Distance  along  flow  passage 

from  the  matrix  inlet 

Z  Reduced  length  (N  x) 

tu  r 


Greek  Letter  Symbols 


P 

A 

e 

A 

A 

J^cmro 

r 


Compactness  (A/V  ) 
tn 

Ratio  of  orifice  diameter  to  pipe 

diameter  (d  /d) 
o 

Difference  or  change  (time, 
temperature,  distance,  etc.) 

Time 

Free  time 
Fluid  viscosity 
Centroid  coordinate 
Density 


ft 

lbm/hr 
lbf/sq  ft 
dimensionless 
Btu/hr 

(ft  lbf)/(lbm  deg  R) 
ft 

deg  F 

dimensionless 

ft/sec 

cu  ft 

lbm 

lbm 

ft 

dimensionless 

sq  ft/cu  ft 
dimensionless 


sec,  hr 
dimensionless 
lbm/hr  ft 
dimensionless 
lbm/eu  ft 
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Subscripts 


atm  Local  atmosphere 

ave  Average 

f  Fluid  (gas,  air) 

i  Initial,  inlet 

m  Matrix,  mean 

o  At  orifice 

s  Solid  (Matrix  material) ,  static 

STD  Standard  (temperature  and  pressure) 

x  Local  conditions 

1  Inlet  conditions  (upstream  of  matrix  and  heaters) 

2  Inlet  conditions  at  matrix  entrance 

3  Exit  conditions  at  matrix  outlet 

Dimensionless  Groupings 

Fanning  friction  factor;  ratio  of  wall  shear 
stress  to  fluid  dynamic  head 

2/3 

Colburn  j-factor  )•  This  factor  plotted  vs. 

Reynolds  Number  defines  the  surface  heat  transfer 
characteristics. 

Longitudinal  heat  conduction  parameter  for  solid 
material  (lc^A^/mLc^) 

Time  parameter  (hA6/Wgc^) 


J 


A 

7r 
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1. 


Introduction. 


The  transient  test  facility  at  the  Naval  Postgraduate  School(NPS) 
has  been  in  operation  for  several  years.  Two  techniques  have  been  used 
at  the  NPS  facility  for  the  determination  of  heat  transfer  data.  One  of 
these  is  the  maximum  slope  technique  developed  by  Locke  [13]  which  uti¬ 
lizes  the  temperature- time  response  curve  of  the  fluid  leaving  the  matrix 
after  a  step  change  in  the  fluid  inlet  temperature.  Howard  [6]  extended 
this  technique  to  include  the  effects  of  longitudinal  conduction.  The 
other  technique  is  the  cyclic  technique  developed  by  Bell  and  Katz  [4], 

In  both  of  these  techniques  the  heat  transfer  parameter  being  sought  is 
N  ,  which  is  a  dimensionless  heat  transfer  parameter  equal  to  the  ratio 
of  the  convective  heat  transfer  rate  from  a  solid  to  the  heating  capacity 
rate  of  an  adjacent  fluid,  i.e., 

where: 


JlA. 

rhsCf 


m  unit  conductance  for  convection  heat 
transfer  (BTU/hr  sq  ft  deg  F) 

*  total  heat  transfer  area  of  solid  (sq  ft) 

*  mass  flow  rate  of  fluid  (lbm/hr) 

*  specific  heat  of  fluid  (Btu/lbm  deg  P) 

It  is  known  that  the  maximum  slope  technique  is  unreliable  at  values 

of  N  less  than  3.5,  due  to  the  large  errors  in  N  associated  with 

tu  °  tu 

errors  in  the  determination  of  maximum  slope.  Furthermore,  it  has  been 
noted  experimentally  that  for  values  of  NCu  near  2.0,  the  temperature 
response  curve  of  a  fluid  displayed  its  maximum  slope  at  approximately 
time  zero  on  the  trace,  and  that  this  value  of  maximum  slope  was  not 
suited  for  determining  N^. 


A 

ft 

m 
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It  was  not  until  the  work  done  by  Kohlmayr  [9,  10],  which  gave  the 
exact  analytical  solution  to  the  single-blow  problem  baaed  on  Hausen's 
[5]  mathematical  model,  that  the  reason  for  these  errors  was  fully  under¬ 
stood.  Later  work  by  Kohlmayr  [11,  12]  demonstrated  ways  in  which  to 
handle  this  problem  as  well  as  the  development  of  a  new  technique  which 
might  be  used  in  place  of  the  maximum  slope  technique  for  values  of  N 

tu 

<5.  Kohlmayr  also  developed  a  means  to  handle  the  single-blow  problem 
when  other  than  s  step  change  in  the  inlet  fluid  temperature  was  made. 

The  purpose  of  this  thesis  was  first  to  determine  experimentally 
the  actual  inlet  temperature  response  of  the  NPS  facility  and  to  use  this 
known  temperature  response  in  the  manner  suggested  by  Kohlmayr.  Addi¬ 
tionally,  a  new  heater  system  was  made  to  try  to  more  closely  approxi¬ 
mate  a  step  change  in  the  fluid  inlet  temperature  so  that  the  physical 
process  and  mathematical  model  might  more  closely  resemble  one  another, 
thereby  improving  the  results  using  the  maximum  slope  technique. 

Another  technique,  the  “time  aero  intercept  technique,"  used  by 
Wheeler  [17]  was  investigated. 
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2.  Summary  of  Theory. 

A.  Background 

The  single-blow  transient  technique  which  is  used  to  determine 
heat  transfer  data  for  a  porous  solid  originally  used  for  its  mathe¬ 
matical  model  Hausen's  [5]  partial  differential  equation  system.  The 
method  involved  comparing  the  recorded  exit  temperature  of  a  fluid  passing 
through  a  porous  solid  which  had  previously  undergone  a  step  change  in  its 
temperature,  with  a  computed  response  curve  based  on  the  solution  to  Hausen's 
equations . 

Locke  [13]  has  shown  that  there  exists  a  unique  relationship 

between  the  maximum  slope  of  the  response  curves  and  the  number  of  heat 

transfer  units,  N  .  However,  Hausen's  model  did  not  include  th  i  effects 
tu 

of  longitudinal  conduction.  Howard  [6],  by  the  use  of  a  finite  differ¬ 
ence  technique  with  the  digital  computer,  included  the  effects  of  longi¬ 
tudinal  conduction. 

It  should  be  noted  here  that  in  all  of  Kohlmayr's  work,  which 
will  be  discussed  later,  the  effects  of  longitudinal  conduction  are  not 
considered . 

B.  Theory 

The  basic  assumptions  in  the  single-blow  problem  are: 

(1)  Properties  of  the  fluid  ere  temperature  independent 

(2)  Fluid  flow  is  steady 

\3)  The  porous  solid  is  homogeneous 

(4)  The  thermal  conductivity  of  both  solid  and  fluid  is 
infinite  in  the  direction  perpendicular  to  flow 

(5)  Thermal  conductivity  of  the  solid  is  zero  in  the 
direction  of  flow 
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Hausen* •  original  differential  equation  system  [5]  was  based  on 
tha  energy  balance  between  a  fluid  passing  through  a  porous  solid  and  the 
solid.  The  system  of  differential  equations  which  resulted  from  this 
development  are  as  follows: 


tf  -  t 
C^ts  ts  =  dr 


in  which  a  "reduced  length"  variable, 
2  = 


and  a  dimensionless  time  parameter, 

■£'=  Q  _  'X  ^ 

K4  m  LjCfXWs^s 

have  been  used. 


In  the  above  equations: 

^5  ■  temperature  of  the  solid  (deg  F) 

“tf  •  temperature  of  the  fluid  (deg  F) 

Q  -  time  (hr) 

“  specific  heat  c<f  the  solid  (Btu/lbm  deg  F) 

Ws  ■  mass  of  the  solid  (Ibm) 

/?L  ■  distance  along  flow  passage  measured  from  inlet  (ft) 

/.  «  total  length  of  solid  (ft) 

*  mass  of  fluid  entrcined  in  colid  (lbm) 

Kohlmayr  (11]  has  modified  Hausen's  original  equations  by  introducing  a 
new  dimensionless  time  variable  called  "free  time," 


u.  =  JDAl  g 

Ws  £s 


W*  </r£- 
ws£t  L 
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(2) 


For  convenience  Kohlmayr  introduced  two  constants. 


»c=  chjk. 


and 


£  = 

'  VJs-Cs. 


which  are  fixed  for  any  particular  experiment.  For  most  practical  pur¬ 
poses,  %  0  and . 

Hausen's  modified  equations  are  then  given  as 

Sz 

where  t*  and  t*  are  the  temperatures  of  the  solid  and  fluid  respec¬ 
tively,  which  have  arbitrarily  been  normalised.  The  above  equations  are 
subject  to  the  boundary  and  initial  conditions: 

«  Z  -  o  tf  =  d(/*}  <5> 

and 

•e/^-o  £$(2,JA)  =  /  (*> 


=  tSM 


where  g( JU.  )  represents  the  normalized,  time-dependent  fluid  tempera¬ 
ture  at  the  inlet  cross-section.  Kohlmayr  solved  equations  (3)  through 
(6)  by  means  of  a  double  Laplace  Transform  [10]. 

The  results  of  Kohlmayr' s  solution  are: 

&(2,a)=  e'2 [g(/A)+  Jz a/6uE, y  g(v)oLv] 


(7) 
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is  an  entire  family  of  functions  of  the  order  K  (K.—  0}  l^. . 


i.e. , 


(j\-\Y 'J)\ 


Evaluated  at 


*  Ll  =  N' 

equation  (7)  becomes 


r»*o 

L 


+ 


L  [/- 

—i lv  yU) (/-  f(-v))o(-v 


(8) 


(9) 


By  differentiation  of  equation  (9)  with  respect  to yU.  , 

+  fQ  w-HLC 

{V*- 


fl  [4  -S/V^a-v;)]  e^^dy, 


(10) 


With  these  results  Kohlmayr  [12]  then  shoved  that  for  other 
than  a  step  change  in  inlet  temperature,  both  the  maximum  slope  and  "free 
time"  at  a  given  value  of  N  may  be  multi-valued.  He  further  demonstrated 
that  by  knowing  what  the  inlet  response  is,  one  can  determine  which  of 
the  values  of  maximum  slope  is  valid  and  also  how  to  determine  what  the 
error  bounds  are  in  the  use  of  some  known  inlet  temperature  change.  One 
of  Kohlmayr's  results  is  that  the  correct  value  of  "relative"  maximum 
slope  occurring  in  the  fluid  exit  temperature  response  will  be  that  value 
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which  occurred  last  in  cases  in  which  Kohlmayr's  technique  is  applicable. 
Furthermore,  the  maximum  slope  method  is  unstable  for  <  2.0,  is 

singular  at  N  ■  2.0,  and  inaccurate  for  Nfcu  between  2.0  and  3.0.  [9,.12l 

In  view  of  the  short-comings  of  the  maximum  slope  technique, 
Kohlmayr  developed  an  indirect  curve -matching  technique  [11].  Some  im¬ 
portant  results  of  this  development  are  summarized  below. 

With  the  solution  to  the  single-blow  problem  known,  it  is  pos¬ 
sible  to  develop  an  indirect  curve -matching  technique  based  on  the  first 
moments  of  the  fluid  transient  response  curves.  This  method,  known  as 
the  centroid  method,  involves  the  reducing  of  both  the  theoretical  tempera¬ 
ture  response  t*(Ntu,yC4)  and  experimental  response  t*^e  (H^,  yM-)  into 
two  different  single-valued  functions  based  on  the  one  parameter  N^. 

To  generalize  the  problem  further,  a  mapping  functional  was  defined, 
based  on  Che  fluid  inlet  temperature  change  and  the  fluid’s  exit  tempera¬ 
ture  response  such  that 


where  the  following  restrictions  are  imposed  on  ; 

(i)  y*  must  be  real,  single -valued,  continuous  with 
respect  to  both  t*(N  )  ant^  )»  an<*  monotone  with  the  parameter 


N 


tu’ 


(2)  (H  )  must  be  monotone  increasing  with  1(g)  * 
,  the  "deviation  from  step." 


cLAk*. 


£  1C 


,  where  K  is  some  measure 


for  the  maximum  permissible  amplification  of  errors. 

(4)  For  any  N£u  and  any  given  deviation  from  step,  (N  ) 
should  be  insensitive  with  respect  to  local  variations  of  g( 


IV 


(5)  The  evaluation  of  (Nfcu)  must  be  simple  and 

straightforward . 

The  functional  chosen  was  the  first  moment  of  the  difference 
between  upstream  and  downstream  fluid  temperatures. 


fj*  —  3+  ~jj 

o  f 


(11) 


Due  to  the  difficulty  involved  with  integrating  equation  (11)  up  to 
values  of  JJ-  ■  Co  ,  a  new  functional  was  chosen  and  defined  as 


ffi  \t*  ( ^  a>-^]  =  ^-+7^ 


(12) 


Some  further  restrictions  imposed  upon  g( yLC )  are: 

(1)  g( JJL  )  must  be  non-negative  monotone  decreasing 

(2)  have  initial  value  g(  O  )  "1  and 

(3)  assume  value  sero  for  all  free  times  which  exceed 
the  maximum  permissible  deviation  from  step:  g( Lk.  ) 

■  0  for  all  Jj_  21  I 

—  max 

Therefore  the  moment  functional  was  defined  in  terms  of  the  ceutroid  of 
the  area  under  this  difference  curve: 


A 


cLjj, 


(13) 
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For  Che  case  of  g( JJ,  )  *  0  (step  change),  the  results  of  equation  (13) 


are  listed  as  Figure  1. 


In  order  to  use  the  above  results  for  any  given  inlet  tempera¬ 
ture  change  g (JJL  ),  Kohlmayr  arrived  at  the  following  empirical  relation¬ 
ship: 


A*J^]  -  A- J^‘0)  +  • 99  A  ±  • 01 


This  relationship  was  obtained  by  systematically  applying  the  centroid 
method  to  many  different  upstream  fluid  temperature  changes. 

Kohlmayr 'a  solution  to  the  single-blow  problem  may  also  be  used 
directly.  Returning  again  to  equations  (9)  and  (1) 

,  u  /,  /  , .  \  /  _  j  r  i  _ .  > 


t:  Wi,  j/0  =  /-  e 


r  i-q-(u-) 


~t~ f  A/ 2  W— y  i 'y)) d  V 


where 


—  m  j~t  &  Ia// 

^  W-  C;  Ws  A;  L_ 

t^(Htu> y6<-  )  can  be  evaluated  at  "free  time,"  *  0,  which  would  corres¬ 

pond  to  the  time  that  an  element  of  fluid  which  has  undergone  a  change  in 
temperature  upstream  arrives  at  the  exit  of  the  solid.  Thus, 


-  ,'A  /  -  3(0) 
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or 


/  ~  K(^0)  =  ^  {'~9(0). 


(16) 


Taking  the  natural  logarithm  of  both  aidea  of  the  above  equation,  yield* 


[/~t^0)]  =  “  C  + 


(17) 


Therefore 


M  = 

tU. 


-£tl  j 


+  JIvl 


I- 


g(0) 


(18) 


Bote  that  fcr  a  step  change  in  the  inlet  temperature,  i,e.,  *(y£4.  ) 
for ^U->0,  that  equatioa  (18)  reduce*  to, 


(19) 


In  order  to  utilize  equation  (19)  directly  it  would  be  necesaary  to  achieve 
a  »tep  change  in  the  inlet  fluid  tesperature,  g(y£/-).  Thi*  would  occur 
at  real  tiae,  Q  *  0.  Then  aeasure  the  temperature  of  that  aaae  element 
of  fluid  which  had  undergone  the  *tep  change  in  temperature,  as  it  left 
the  porous  solid.  This  would  occur  at  "free  tiae,"  LA-  *  0. 
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Because  of  the  finite  response  tine  of  the  different  components 
of  the  experimental  equipment,  i.e.„  heaters,  thermocouples  and  recorder, 
it  is  not  possible  to  aeet  precisely  cither  of  the  above  requirements. 

In  viev  of  the  above  limitations  it  is  necessary  either  to  make  some  ap¬ 
proximations  in  the  interpretation  of  the  temperature  response  data,  t£, 
or  else  devise  a  means  by  which  the  actual  response  data  might  be  treated 
directly.  Some  of  the  difficulties  involved  with  the  second  alternative 
are:  Refer  again  to  equation  (9) 


+ ■/fczfiQti'-M cLv\ 

which  nay  be  stated  as 

l  I 

+ //£  e'1-’ 

O 

and  the  cxpresaion  for  the  inlet  fluid  temperature, 

=  , 


where  0M  is  the  equal  to  the  heater  time  constant  which  was  derived  in 
Appendix  C. 

One  approach  to  using  equation  (9)  would  be  to  seasure  t*  at 
some  value  of  "free  tiae'yU.  >  0  such  that  g (yU.  )  =  0  for  which 
would  be, 


At  = 

/  +  ^  <r£j.  v-yu))  e~Mi^V/U)0-gi  v))d  y 


The  difficulty  here  is,  that  in  order  to  evaluate 


23 


H  mu»t  first  be  known.  This  would  result  in  a  direct  curve  raatchina 
tu 

technique.  A  second  approach  would  be  to  measure  t£  at  some  time  0 

such  that  /M- 

Q(p.)  =  /il 

Jo 

again  Ntu  must  first  be  known.  In  both  the  above  methods  if  yi/L  is  small 
the  transient  response  of  the  sensing  and  recording  equipment  will  still 
present  problems. 

Because  of  the  difficulties  encountered  with  the  above  two 
method,  *:he  first  alternative,  that  of  making  some  approximations  in  the 
interpretation  of  the  measured  response  data,  t£,  will  be  considered. 

Since  the  response  of  the  thermocouples  is  taster  than  that  of  the  .003" 
diameter  heaters,  and  approximately  equal  to  that  of  the  .001"  heaters, 
see  Appendix  C,  it  might  be  reasonable  to  assume  that  the  initial  reapon  e 
recorded  during  a  particular  run  is  that  response  due  to  the  heaters,  see 
Figure  2B.  Therefore,  it  is  possible  to  extrapolate  back  to  "free  time" 
/''•  ■  0,  from  the  tii.e  on  the  response  curve  where  the  transient  response 
of  the  heaters  has  died  out  and  treat  the  intersection  of  the  extrapolated 
curve  and  ,vfree  time,"^/-  *  0  ns  an.  actual  step  change  in  temperature  The 
above  technique  is  referred  to  as  the  "zero  intercept"  technique. 
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3.  Experimental  Technique. 

The  existing  apparatus  at  the  NPS  facility,  aet  Figureo  3  and  4, 
has  been  desigued  to  conform  to  the  idealisations  required  by  Howard  [5] 
for  the  use  of  his  conduction  parameter  in  the  maximum  slope  technique. 
Howard's  conduction  parameter,  ^  ,  is  defined  as 

^  =  -4s  fts 

Yr\  LjC+ 

where : 

-4$  *  thermal  conductivity  of  the  solid  (Btu/hr  aq  ft  deg  F/ft) 

As  «  matrix  solid  cross-sectional  area  available  for  thermal 

conduction  (sq  ft) 

Howard's  idealisations  are  as  follows: 

(1)  Hie  fluid  flow  in  the  matrix  is  both  steady  and  uniform 
in  velocity  and  temperature  at  any  cross  section 

(2)  The  matrix  thermal  conductivity  is  finite  in  the  direction 
parallel  to  fluid  flow  and  infinite  in  the  direction  normal  to  flow 

(3)  The  matrix  thermal  conductivity  is  large  in  comparison  to 
that  of  the  contained  fluid 

(4)  The  thermal  properties  of  the  fluid  and  matrix  are  conatant 
and  uniform 

(5)  The  convective  heat  transfer  coefficient  is  some  suitable 
average  and  remains  constant 

(6)  A  step  change  in  the  temperature  of  the  inlet  fluid  is  im¬ 
posed  at  real  time  equal  to  zero. 

These  idealizations  result  in  making  the  fluid  flow  one  dimensional. 
They  in  no  way  conflict  with  the  restrictions  imposed  by  Kohlmayr  [11], 
The  requirement  of  uniform  velocity  and  temperature  profiles  (1) 
is  met  by  a  specially  designed  entrance  nozzle,  flow  straightening 
screens,  and  an  even  distribution  o£  heater  wires  across  the  channel. 
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Piereall  [15],  using  the  equipment,  verified  that  the  velocity  and  tempera 
ture  profiles  were,  indeed,  uniform  within  reasonable  limits.  The  im¬ 
portance  of  uniform  velocity  and  temperature  profiles  was  demonstrated 
by  Wheeler  [17], 

The  ■  .  temperature  change  across  the  matrix  of  about  20®F  gives 

constant  thermal  properties  of  the  fluid  and  matrix  and  a  constant  heat 
transfer  coefficient  as  required  by  idealisations  (5)  and  (6). 

The  rastriction  that  the  temperature  of  the  inlet  fluid  be  subjected 
to  a  step  change  will  be  discussed  later. 

The  data  roquired  for  the  computation  of  the  various  heat  transfer. 


fluid  flow  coefficients,  and  dimensionless  parameters  is  as  follows: 
fariA  *  atmospheric  pressure  (wo  Hg) 

■  orifice  static  pressure  (inches  IIjO) 

•  pressure  drop  acrcss  orifice  (inches  H^O) 

^  -  pressure  drop  across  matrix  (inches  H20) 

py  ■  static  inlet  pressure  at  entrance  to  heaters 
(Inches  I^C) 

"to  »  temperature  of  fluid  at  orifice  (millivolts) 
do  ■  orifice  diameter 


P 

cs 


*  ratio  of  orifice  diameter  to  pipe  diameter 

*  chart  spe'ad  (sec/inch) 


ts-t, 

t3-tl 


-  downstream  response 

-  downstream  -  upstream  response 


Pressures  are  measured  with  either  a  draft  gage  or  a  water  manometer, 


depending  upon  the  orifice-flow  rate  combination,  except  for  atmospheric 


pressure  which  is  measured  by  a  mercury  barometer. 

The  temperature  responses  recorded  are  t^-t^,  the  difference  between 
the  inlet  fluid  temperature  and  the  fluid  temperature  at  the  matrix  exit 


(used  'n  the  maximum  slope  and  zero  intercept  techniques)  and 
the  difference  between  the  fluid  exit  temperature  and  the  matrix  inlet 
fluid  temperature  (used  in  the  centroid  technique).  For  a  more  complete 
description  of  equipment,  the.  reader  is  referred  to  Appendix  A.  See 
Figure  5  for  the  position  of  temperature  and  pressure  measurements. 

A  test  run  is  accomplished  by  predetermining  the  necessary  pressure 
drop  across  the  orifice  to  achieve  a  desired  flow  rate  and  to  determine 
the  number  of.  heaters  necessary  to  achieve  a  20°F  temperature  rise.  Air 
is  drawn  through  the  apparatus  and  is  controlled  at  the  entrance  to  the 
turbocompressor.  When  the  desired  flow  rate  is  achieved  the  heaters  are 
then  energized  and  the  heated  air  and  test  core  are  allowed  to  reach  a 
steady  state  temperature,  at  which  time  pressure  measurements  are  re¬ 
corded.  Power  to  the  heaters  is  then  secured  and  recordings  of  t^-^ 
and  t3~t2,  both  as  functions  of  time  are  recorded  on  separate  channels 
of  the  Brush  recorder  (see  Figures  2A  and  B).  The  temperature  at  the 
orifice  (t^)  is  measured  before  and  after  a  run  to  insure  the  same  ambient 
air  temperature  for  the  run. 

After  completion  of  the  desired  runs,  the  values  needed  to  compute 
slope  are  taken  from  th)  recorded  traces  of  t^-C^.  The  maximum  value  of 
slope  is  obtained  visually  with  the  aid  of  a  straight  edge.  This  informa¬ 
tion  is  included  on  the  data  sheet.  The  data  sheet  layout  conforms  to 
the  data  input  section  of  a  digital  computer  program  [18],  which  reduces 
the  data  used  ir  the  maximum  slope  technique  and  calculates  the  parameter 
C?C  used  in  the  centroid  technique.  For  complete  details  on  data  re¬ 
duction  the  reader  is  referred  to  Appendix  B  at  this  time. 
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4.  Description  of  Test  Matrices. 

The  two  cores  used  In  this  experiment  were  both  cores  which  hsd 
been  tested  previously  end  for  which  there  was  a  good  deal  of  data  avails 
able. 

One  matrix,  the  Cercor  T20-38,  previously  tested  by  Howard  [7],  ia 
a  ceramic  type-core  having  a  low  thermal  conductivity,  and  was  used  to 
eliminate  any  errors  that  might  be  Introduced  by  longitudinal  conduction 
which  might  have  invalidated  results  using  the  centroid  and  aero  inter¬ 
cept  techniques.  The  other  core  tested,  which  is  a  stainless  steel  plate- 
fin  type  matrix.  Solar  4,  had  been  tested  several  times  [15],  [16],  [18] 
and  was  used  so  that  comparisons  could  be  made  when  longitudinal  conduc¬ 
tion  was  a  factor. 

Further  information  on  core  geometries  and  properties  is  shown  in 
Figures  6  and  7. 
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5.  Presentation  of  Results. 

For  each  matrix  tested,  the  heat  transfer  and  flow  friction  character¬ 
istics  have  been  computed.  The  computed  results  are  shown  tabulated  in 
Tables  I-II.  The  values  of  Nfcu  have  been  plotted  for  the  three  methods 
used  in  this  experiment.  Figures  8A  and  9  compare  values  of  M£u  based  on 
the  maximum  slope  and  centroidal  technique  for  cores  T20-2o  and  Solar  4. 
Figuren  8B  and  10  compare  values  of  N  based  on  the  maximum  slope  and  zero 
intercept  technique  for  the  same  two  cores.  Figure  11  is  a  comparison  of 
N  values  computed  by  the  maximum  slope  technique  for  the  two  heater 
systems.  Core  Solar  4  was  the  only  core  tested  with  the  new  .001  inch 
diameter  wire  heaters. 

The  heater  response  curves  are  shown  in  Figure  12.  See  Appendix  C 
for  complete  details  used  in  determining  these  curves. 
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6.  Discussion  of  Results 

In  order  to  insure  that  the  technique  used  by  the  author  in  this 
experiment  was  correct,  the  results  of  the  single-blow  test  of  this  experi¬ 
ment  were  compared  with  the  results  of  previous  tests  of  the  same  cores. 

Considerable  difficulty  was  experienced  in  the  use  of  the  new  (.001 
inch  diameter)  heater  system.  The  heaters  burned  out  frequently,  thereby 
limiting  the  range  of  Reynolds  Numbers  over  which  the  core  was  tested  and 
also  the  number  of  cores  actually  tested.  The  exact  reason  for  the  burn¬ 
out  is  not  known.  It  was  initially  thought  to  be  due  to  surges  in  cur¬ 
rent  caused  by  the  use  of  a  resistor  bank  in  the  first  experiments.  How¬ 
ever,  when  the  resistor  bank  was  removed  and  voltage  to  the  heaters  control¬ 
led  by  the  motor  generator  set  rheostat,  the  system  failed  again.  In  all, 
there  were  three  different  failures,  all  occurring  at  different  flow  rates 
and  with  a  different  number  of  heaters  in  ure  at  the  time  of  failure. 

An  examination  of  Figure  11  shows  that  in  the  high  Reynolds  number 
ranges  there  is  little  difference  in  the  value  of  N  based  on  the  maximum 
slope  data  for  the  two  different  heater  systems.  Furthermore,  this  is  in 
a  range  where  the  errors  in  Nfcu  due  to  errors  in  the  maximum  slope  techni¬ 
que  are  greatest.  Without  actually  applying  Kohlmeyr's  equations  in  the 
"extension  of  the  maximum  slope"  technique  [17],  it  is  impossible  to  get 
a  quantitative  value  for  the  error  in  N  due  to  deviations  from  the  step 
change.  However,  by  the  use  of  his  curves  based  on  the  "deviation  from 
step"  (previously  defined  as  1(g)  and  equal  to  the  area  under  the  inlet- 
temperature  response  curve)  which  are  presented  here  as  Figure  13,  it  can 
be  seen  that  for  I  <<.100  the  error  in  slope,  Am  which  is  the  difference 
in  the  maximum  slope  due  to  a  step  change  in  the  inlet  fluid  minus  the 
maximum  slope  due  to  a  non-step  change  in  the  inlet  fluid,  is  less  than 
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.01.  It  must  be  pointed  out  that  Figure  13,  is  for  an  upstream  tempera- 

2 

ture  change  which  is  quadratic,  i.e.,  g (jj.  )  *  (1-yW/A)  .  The  curves  of 

Am  versus  N  for  the  NPS  facility  would  be  slightly  different  from 
t  u 

those  presented  In  Figure  13,  for  the  g (JJL)' z  associated  with  the  NPS 
facility's  heate:  systems  are  exponential.  However,  for  low  values  of 
I(I<  .  100)  the  error  in  maximum  slope  Am  would  be  small  regardless  of 
the  exact  form  of  the  inlet  temperature  responsi.  Since  the  maximum  1(g) 
resulting  from  either  heater  system  is  less  than  .03  It  would  not  be 
expected  fo*  one  to  detect  any  difference  of  slope  and  consequently  N 
resulting  from  the  use  of  either  system. 

Again  referring  to  Kohlmayr's  curves.  Figure  13,  deviations  from 
step  changes  In  temperatures  would  cause  large  errors  at  lower  flow  rates, 
l.e,,  Ntu>5.  However,  since  the  1(g)  of  the  NPS  facility  decrease* 
with  decreasing  flow  rates  (see  Appendix  C)  it  could  be  assumed  that  the 
errors  in  maximum  slope  due  to  1(g)  would  be  small.  Furthermore,  since 
the  longitudinal  conduction  is  greater  at  low  flow  rates  its*  effects 
cannot  be  neglected  as  they  have  in  Kohlmayr's  assumptions. 

In  the  Investigation  of  the  centroid  and  the  zero  intercept  techni¬ 
ques  the  range  of  flow  rates  of  both  cores  tested  was  the  same  snd  ranged 
from  m  ~  250  to  ff)  v>  950  Ibm/hr.  The  reason  for  950  lbm/hr  being 

the  maximum  value  tested  was  that  at  this  flow  rate  all  heaters  were  in 
use  and  in  order  to  go  to  higher  flow  rates,  an  increase  in  voltage  would 
have  been  necessary  ro  achieve  the  same  20  deg  F.  temperature  rise.  This 
would  ha  e  changed  the  beater  time  constant,  which  is  a  function  of  both 
flow  rate  and  the  number  of  heaters  in  use. 

Investigation  of  the  curves  of  Figure  8A  and  9  reveals  that  for  both 
cores  teoted  there  is  very  close  agreement  between  the  results  predicted 
by  the  maximum  slope  technique,  shown  as  the  dotted  line,  and  the  results 
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of  the  centroid  technique  for  values  of  N  between  3.5  and  5.  In  fact 
there  is  even  some  overlapping  of  points  in  Figure  8A.  For  values  of 
Ntu  below  3.5,  the  values  of  N  based  on  the  maximum  slope  technique 
begin  to  decrease  from  the  predicted  values  whereas  the  values  determined 
using  the  centroid  technique  followed  closely  the  predicted  values  for 
the  case  of  Solar  4  and  increased  slightly  from  the  predicted  values  for 
core  T20-38. 

It  can  be  seen  from  Figure  8B,  which  compares  the  values  of  H 
computed  by  the  maximum  slope  and  zero  intercept  method,  for  core  T20-38, 
that  the  flow  rates  used  were  not  high  enough  for  the  zero  intercept 
method  to  be  used  to  its  best  advantage.  Only  those  runs  made  at  flow 
rates  greater  than  600  lbrn/hr  had  a  clearly  measurable  sero  intercept. 

Figure  10  again  compared  maximum  slope  N^'s  with  those  evaluated 
by  the  aero  intercept  technique  in  this  case  for  core  Solar  4.  This  was 
done  for  both  sets  of  heaters.  It  is  important  to  note  the  close  correla¬ 
tion  between  the  values  of  Ncu  for  the  two  heater  systems.  This  will  be 
covered  more  fully  in  Section  7. 
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7.  Experimental  Uncertainties. 

Various  idealizations  and  boundary  conditions  have  been  imposed  for 
the  mathematical  model  of  the  physical  experiment.  Due  to  the  fact  that 
these  idealizations  and  boundary  conditions  have  not  been  precisely  met, 
certain  errors  have  been  introduced.  These  errors  are  difficult  to  as¬ 
sign  a  numerical  value  to,  and  with  the  exception  of  the  deviation  from 
step  temperature  change  previously  mentioned  and  the  effects  of  longi¬ 
tudinal  conduction,  will  not  be  discussed. 

The  experimental  errors  associated  with  the  maximum  slope  technique 
used  at  the  NFS  facility  have  been  covered  in  considerable  detail,  [2j[3] 
f 15 3 ;  therefore  only  those  errors  introduced  in  this  particular  experi¬ 
ment  will  be  discussed. 

In  using  Kohloayr's  centroidal  technique  there  are  several  sources 
of  possible  error.  First,  Kohlmayr,  in  using  Hausen's  mathematical  model 
has  neglected  longitudinal  conduction  in  the  solid.  Howard  [6]  has  calcu 
lated  the  errors  in  N  associated  with  longitudinal  conduction  in  using 
the  maximum  slope  technique,  but  no  such  information  is  available  for  use 
with  the  centroidal  technique,  for  the  analytical  solution  to  the  single¬ 
blow  problem,  in  this  case,  does  not  include  the  effects  of  longitudinal 
conduction.  Therefore,  one  can  only  get  a  qualitative  idea  as  to  the  ef¬ 
fects  of  these  errors. 

Fi'om  Howard's  curves  of  N  versus  maximum  slope.  Figures  14  and  15, 
for  a  given  value  of  N  the  slope  of  a  temperature  response  curve  is 
different  for  different  values  of  the  conduction  parameter  A  .  There¬ 
fore,  the  response  curves  and  consequently  the  centroid  of  the  response 
curves  would  be  different  for  different  values  of  A  .  The  difference  in 

maximum  slope  due  to  A  decreases  as  N  decreases.  For  example  at  N 

tu  r  tu 

30  a  difference  in  A  of  approximately  .1  causes  a  90%  difference  in 


33 


maximum  slope,  whereas  for  Ntu<5  the  maximum  error  in  maximum  slope 

caused  by  *  .1  difference  in  A  is  approximately  9%.  It  would  seem 

reasonable  Co  assume  that  provided  Kohlmayr's  results  were  applied  only 

in  the  low  N^u  range  (Ntu<5)  and  longitudinal  conduction  was  small  (i.e., 

A  <  .1),  conduction,  errors  could  be  neglected.  For  this  experiment 

the  maximum  value  of  A  for  either  core  at  values  of  N,  <5  was  less 

tu 

than  .02?. 

Another  source  of  error  is  in  the  determination  of  the  centroid  of 
the  area  under  the  response  curve.  Two  methods  were  used  for  this  experi¬ 
ment,  primarily  as  checks  on  each  other.  In  both  methods  the  critical 
point  is  the  determination  of  where  the  cut-off  is  for  a  particular 
curve.  The  cut-off  point  is  that  point  on  the  response  curve  which  is  one- 
tenth  of  the  maximum  value.  This  point,  defined  as  occurs  at  a 

point  where  the  slope  of  the  curve  is  slowly  approaching  »ero,  and  an 
error  of  one  millimeter  in  the  ordinate  can  cause  an  error  of  several 
millimeters  in  the  abscissa.  For  example,  in  the  data  taken  with  core 
Solar  4  when  Ntu  »  1.41,  at  (t^-t^)  *  kam,  0  *  13.0  sec.,  and  for  (t^-t^) 

■  5,  0  *13.75  aec.  This  was  for  a  typical  trace  in  which  (t^-t2)max  * 

50.5mm.  A  sketch  of  the  response  curve  (see  Figure  16)  will  help  to 
clarify  the  problem. 

The  actual  error  in  y^cENTD  causec*  ^y  an  error  determining 
will  depend  on  the  value  of  Ntu.  A  typical  example  would  be  for  Nfcu  - 
3.5.  When  computed  by  the  manual  technique  (see  Appendix  B  for  details), 
it  was  found  that  a  +  1  tan  difference  in  ordinate  (which  is  ^  2%  and 
would  be  a  maximum)  caused  an  errer  of  2T»  ^ny^0gjjTD'  8*10u*<* 
pointed  out  that  these  values  are  approximate,  for  there  is  some  uncertain¬ 
ty  aasociated  with  finding  the  intersection  of  the  three  centers  of  mass 
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lines,  shown  in  Figure  16  and  discussed  in  Appendix  B.  This  uncertainty 
usually  was  less  than  +  27..  The  same  type  of  error  associated  with  find¬ 
ing  the  cutoff  point  is  also  found  in  the  computer  technique  of  determining 
the  centroid.  In  this  method  the  computer  solution  required  data  points 
from  the  response  curve  at  fixed  intervals,  which  for  this  experiment  were 
5mm  intervals.  Therefore,  if  the  value  of occurred  at  a  point  mid¬ 
way  in  any  one  5rm  increment,  there  could  be  as  much  as  a  2%  mm  error  in 
the  cut-off.  For  an  average  130  to  140  can  trace  this  is  less  than  2%. 

For  comparative  purposes,  data  runs  were  *.<ide  using  both  techniques  and 
the  maximum  deviation  in  centroid  between  the  two  methods  was  found  to  be 
4%,  but  the  majority  of  deviations  v„s  less  than  11. 

Another  difficulty  discovered  in  analyzing  the  temperature  trace  was 
caused  by  the  sensitivity  of  the  thermocouples.  During  a  particular  run 
a  change  in  the  steady  flow  caused  by  a  sudden  draft  in  the  laboratory 
caused  a  temperature  deviation  on  the  trace.  This  required  some  visual 
smoothing  on  the  part  of  the  author  to  extrapolate  the  actual  response. 

This  might  cause  errors  when  using  an  automated  data  reduction  process, 
unless  it  included  a  good  smoothing  technique. 

An  additional  problem  was  encountered  in  the  process  of  data  reduc¬ 
tion.  One  of  the  assumptions  made  in  using  the  centroid  technique  *  that 
at  time  zero  the  temperature  both  upstream  and  downstream  from  the  matrix 
is  constant.  It  was  found  that  this  was  not  the  case  for  this  experiment, 
the  reason  being  that  there  was  some  heating  of  the  thermocouples  located 
adjacent  to  the  heaters  due  to  radiation  from  the  heaters.  This  meant 
that  the  zero  reference  point  of  the  trace  of  e*tablishea  before 

and  after  the  run  (heaters  off)  was  not  the  same  as  that  just  prior  to 
time  zero  when  the  heaters  were  deenergized  (see  Fig.  16).  Therefore, 
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where  the  response  curve  8hou)d  have  resembled  curve  a  of  Fig.  16,  it 
took  the  form  of  curve  b.  To  compensate  for  this,  point  c  was  taken  as 
time  zero  and  curve  a  used  to  determine  the  centroid  It  should  be  pointed 
out  here  that  these  figures  are  greatly  exaggerated  here  for  clarity  and 
the  error  caused  by  this  difference  is  felt  to  be  negligible.  Also,  this 
difference  was  most  evident  at  low  flow  rates,  f or  the  radiation  effect 
was  caused  primarily  by  che  heaters  nearest  the  thermocouples.  As  the  flow 
rates  Increased  and  more  heaters  were  energized  the  percentage  of  heating 
due  to  radiation  became  less. 

Sc  far  no  mention  has  been  made  of  the  error  in  N  caused  by  an 
error  in  /*■  .  Kohlmayr  [11]  has  conducted  a  linear  error  analysis 

based  on  the  approximation 

A  Mtou  *  ciJ^~  AAcea ,to 
c^/lce/^to 

and  introduced  a  relative  error  amplification  factor 

«  £ 

d/t-CEWTD  U*- 

or 

~  K  ^centd 

With  empirical  results  he  has  tabulated  values  of  K  versus  ancI 


K  versus  IT  . 

tu 

As  an  example,  using  Kohlmayr' s  curves  when/^cg^xD  ia  ®ea8ured  t0° 
high  by  27,  and  when  K  -  -5.0,  tnen  N  as  obtained  by  this  method  would 
he  10%  too  low. 

Associateo  with  the  zero  intercept  technique  are  two  distinct  types 
of  error.  One  reeults  from  the  initial  assumptions  on  which  this  technique 
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is  based.  These  assumptions  were  discussed  in  Section  2.  If  errors  are 
present  in  the  results,  they  may  be  due  to  errors  in  the  original  assump¬ 
tions.  It  is  not  possible  to  get,  at  this  tine,  a  quantitative  value  for 
this  type  of  error  or  in  fact  to  even  verify  that  the  technique  used  here 
is  valid.  The  only  fictor  which  would  substantiate  the  technique's  vali¬ 
dity  is  that  for  both  heater  systems  used  in  the  experiment  the  results 

in  terms  of  N  were  nearly  identical.  The  second  type  of  error  is  that 
tu 

associated  with  the  physical  data  reduction,  primarily  the  correct  extra¬ 
polation  of  the  response  curve  back  to  tine  zero,  sef  Figure  17.  This 
may  be  seen  in  the  following  development: 

Mtu_  -  -  JU(l  -  t*) 


y^±  =  a  if/) 


=  )-t*  ±  At? 


e 


tu. 


=  /±  AM 


/-£* 


or 


+  =  a±  a ti_  \ 

'-*? 

so  that  if  1-t*  *  .2  corresponding  to  a  value  of  N  *•  1.609,  and 
f  tu 


a 

-+.oz 

or  2% 

> 

II 

•  ■+.  OZOZ 

or  1,25% 

- .  oz. 

or  -27. 
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0/98 


or  -1.2%. 


8. 


Cone lusions . 


1.  The  present  design  of  the  new  .001  inch  diameter  heaters  was 
not  found  to  be  satisfactory  due  to  th-ir  frequent  and  sometimes  un¬ 
explained  failures.  No  improvement  of  the  results  of  the  maximum  slope 
and  zero  intercept  tests  was  obtained  using  the  smaller  diameter  wire 
heaters.  No  results  were  obtained  with  the  new  heaters  for  evaluation 
with  Che  centroid  technique,  but  it  should  be  pointed  out  that  for  the 
centroid  technique  it  is  not  necessary  to  have  a  step  change  in  the  in¬ 
let  temperature  so  long  as  the  actual  change  is  known. 

2.  It  would  be  extremely  difficult  to  try  and  compare  the  three 

techniques  used  in  this  experiment  on  the  basis  of  test  results  ac¬ 
curacy.  The  best  that  one  can  do  at  present  is  to  compare  the  various 
experimental  results  with  the  predicted  values  based  on  the  maximum 
slope  technique  evaluated  at  high  values  of  N  (Nfcu  >  5.0).  In  this 
respect  the  maximum  slope  technique  is  unreliable  for  Ntu<C  3.5.  The 
centroid  technique  appears  to  have  good  results  for  5  to  as  low 

as  .75  which  was  the  minimum  value  tested  in  this  experiment.  For  both 
the  cores  tested  using  the  centroid  technique  the  results  were  either  the 
same  as  those  predicted  by  the  maximum  slope  technique  or  somewhat  higher. 
The  zero  intercept  method  appears  to  be  impractical  for  values  of  Ntu^>2.5 
due  to  the  larger  differences  associated  with  the  logarithm  of  very  small 
numbers.  This  might  be  better  understood  by  looking  at  the  governing 
equation  in  the  zero  intercept  technique 

Mtu.  =  —  Jho 

and  noting  that  as  the  argument  of  the  logarithmic  term  gets  smaller  the 
natural  logarithm  itself  approaches  negative  infinity.  Below  the  value  of 
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Nt,u<2.5  the  results  were  higher  than  predicted  but  were  extremely  con¬ 
sistent,  that  is,  for  different  runs  conducted  at  the  same  conditions 
the  results  were  nearly  identical. 

From  the  standpoint  of  ease  of  evaluating,  the  zero  intercept  techni 
que  is  very  fast  and  simple  to  apply.  Furthermore,  in  the  method  used  in 
this  experiment  the  zero  intercept  technique  appears  to  be  insensitive 
to  slight  deviations  from  a  step  temperature  change.  On  the  other  hand 
the  results  of  the  centroid  technique  rel^  entirely  upon  how  well  the  so 
called  "deviation  from  step"  change  of  the  fluid's  upstream  temperature 
is  known. 

The  centroid  technique  is  much  more  ‘;ime  consuming  to  use  and  in 
view  of  the  fact  that  its  results  are  based  on  a  considerable  portion  of 
the  fluid  temperature  response  trace  it  is  more  subject  to  errors  caused 
by  sudden  fluctuations  in  the  ambient  conditions,  for  in  the  zero  inter¬ 
cept  technique  only  the  very  first  part  of  the  trace  is  of  interest. 

Based  upon  the  above  conclusions  the  following  table  is  recommended  as 
a  guide  in  determining  which  of  the  aforementioned  techniques  should  be 
used  for  a  given  test  range. 


AO 


9. 


Recommendations  for  Further  Study. 


It  is  recommended  that  an  attempt  be  made  to  determine  both  experi¬ 
mentally  and  analytically  the  effects  of  deviations  from  experimental 
assumptions  such  as  constant  fluid  properties  and  the  convective  heat 
transfer  coefficient  on  the  results  using  the  different  techniques 
including  the  cyclic  testing  technique. 


41 


10.  Bibliography. 


1.  American  Society  of  Mechanical  Engineers.,  Power  Teat  Code 
Supplements,  Instruments  and  Apparatus,  PTC-19.5;  4,  Flow 
Measurements,  Chap.  4,  New  York:  ASME  1959. 

2.  Ball,  Stuart  F.,  "Experimental  Determination  of  Heat  Transfer 
and  Flow  Friction  Characteristics  for  Several  Plate  Fin  Type 
Compact  Heat  Exchanger  Surfaces,"  Master's  Degree  Thesis, 

Naval  Postgra  ‘  late  School,  Monterey,  California,  1966. 

3.  Bannon,  John  M. ,  "An  Experimental  Determination  of  Heat  Transfer 
and  Flow  Friction  Characteristics  of  Perforated  Material  For 
Compact  Heat  Exchanger  Surfaces,"  Master's  Degree  Thesis, 

Naval  Postgraduate  School,  Monterey,  California,  1964. 

Bell,  J.  C.  and  Katz,  E.  F.,  "A  Method  for  Measuring  Surface 
Heat  Transfer  Using  Cyclic  Temperature  Variations,"  Heat  Trans¬ 
fer  and  Fluid  Mechanics  Institute,  Berkeley,  California,  1949. 

5.  Hausen,  H. ,  Z.  Math,  Mech.  9,  page  173,  1929. 

6.  Howard,  C.  ?. ,  "The  Single  Blow  Problem  Including  The  Effects 
of  Longitudinal  Conduction,"  ASME  paper  number  64-GTP-ll,  1964. 

7.  Howard,  C.  P. ,  "Heat-Transfer  and  Flow  Friction  Characteristics 
of  Skewed -Pas sage  and  Glass-Ceramic  Heat-Transfer  Surfaces," 

TR  No.  59,  Department  of  Mechanical  Engineering,  Stanford 
University,  Stanford,  California,  October  1963. 

8.  Kays,  W.  M. ,  and  London,  A.  L. ,  Compact  Heat  Exchangers, 

Second  Edition,  McGraw-Hill  Book  Co.,  Inc.,  New  York,  1964. 

9.  Kohlmayr,  G.  F.,  "Exact  Maximum  Slopes  for  Transient  Matrix 
Heat  Transfer  Testing,"  Pratt  &  Whitney  Aircraft,  Report  FWA- 
2641. 

10.  Kohlmayr,  G.  F.,  "Analytical  Solution  of  the  Single-Blow 
Problem  by  a  Double  La  Place  Transform  Method,"  Pratt  &  Whitney 
Aircraft,  Report  PWA-3158,  July  1967. 

11.  Kohlmayr,  G.  F.,  "An  Indirect  Curve  Matching  Method  for  Trans¬ 
ient  Matrix  Heat-Transfer  Testing  in  the  Low  N  Range,"  Pratt 
&  WMtney  Aircraft,  Report  PWA-3207,  September  u1967. 

12.  Kohlmayr,  G.  F.,  "Extension  of  the  Maximum  Slope  Method  to 
Arbitrary  Upstream  Fluid  Temperature  Changes,"  Pratt  &  Whitney 
Aircraft,  Report  WA-3228. 

13.  Locke,  G.  L.,  "Heat  Transfer  and  Flow  Friction  Characteristics 
of  Po.ous  Solids,"  TR.  No.  10,  Department  of  Mechanical  Engine¬ 
ering,  Stanford  University,  Stanford,  California,  1  June  1950. 


42 


14.  Murdock,  J.  W.,  "Tables  for  Interpolation  and  Extrapolation 
of  ASME  Coefficients  for  Square-Edged  Concentric  Orifices," 

ASME  Paper  No.  64  -WA/FM-6. 

15.  Piersall,  C.  H.,  "Experimental  Evaluation  of  Several  High 
Performance  Surfaces  for  Compact  Heat  Exchangers,"  Master's 
Degree  Thesis,  Naval  Postgraduate  School,  Monterey,  California, 
1965. 

16.  Traister,  R.  E.,  "Experimental  Evaluation  of  Heat  Transfer 
and  Flow  Friction  Characteristics  of  Several  Compact  Heat  Ex¬ 
changer  Surfaces  Utilizing  the  _.nglc  Blow  and  Cyclic  Methods," 
Master'd  Degree  Thesis,  Naval  Postgraduate  School,  Monterey, 
California,  1967. 

17.  Wheeler,  A.  J. ,  "Single-Blow  Transient  Testing  of  Matrix- 
Type  Heat  Exchanger  Surfaces  At  Low  Values  of  Ntu,"  TR.No. 

68,  Department  of  Mechanical  Engineering,  Stanford  University, 
Stanford,  California,  May  1968. 

18.  Trost,  H.  J. ,  Jr.,  "The  Use  of  The  Analog  Computer  With  the 
Single  Blow  Transient  Testing  Technique  For  Compact  Heat  Exchanger 
Surfaces,"  Master's  Degree  Thesis,  Naval  Postgraduate  School, 
Monterey,  California,  September  1967. 


43 


APPENDIX  A 


Description  of  Equipment 

Air  Supply. 

The  working  fluid  for  the  experiment  is  air  drawn  through  the  test 
equipment  by  a  30  HP,  multi-stage  Spencer  Turbo-Compressor,  which  is 
rated  at  550  cfm  operating  on  a  220  V  a.c.  power  supply  see  Figure  3. 

Flow  Measuring  System. 

Flow  measurement  was  accomplished  with  an  ASME  standard  orifice 
section.  Pressure  taps  were  located  d  and  d/2  diameters  upstream  and 
downstream  respectively  from  the  orifice.  Thin  concentric  orifices  with 
throat  diameters  of  .775,  1.232,  1.540,  and  2.310  inches  were  used. 

Heater  System. 

Two  different  heater  systems  were  used.  One  system  utilised  a  .0031 
inch  diameter  nichrome  wire  as  the  heating  element,  and  the  other  .001 
inch  nichrome  wire.  Both  systems  were  designed  to  give  the  same  nominal 
resist* uce  for  a  given  number  of  heater  switches. 

The  .0031  inch  system, (Figure  18A) ,  consisted  of  14  separate  bakelite 
frames,  each  wound  with  two  parallel-connected  heater  elements.  Each 
pair  of  heaters  was  controlled  by  an  individual  switch,  which  in  turn  was 
wired  in  parallel  with  the  other  heater  switches.  All  were  controlled 
by  one  master  switch. 

The  .0010  inch  system  (Figure  18B)  consisted  of  the  same  number  of 
bakelite  frames,  with  each  frame  containing  six  heater  elements  connected 
in  parallel.  The  number  of  heaters  and  voltage  could  be  varied  to  achieve 
approximately  a  20*F.  temperature  change  for  any  given  flow  rate. 

Matrix  Holder  and  Test  Section. 

The  matrix  holder  and  test  section  (Figure  4)  *ere  both  constructed 
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of  polyethylene  plastic.  The  matrix  holder  is  a  drawer  which  slides  into 
the  test  section.  A  removable  frame  on  the  side  of  the  holder  allows  for 
the  insertion  of  the  matrix.  The  flow  channel  is  3-1/16  inches  by  3-1/16 
inches,  and  can  hold  matrices  up  to  3  inches  long.  The  matrix  is  held 
snugly  in  place  by  styrofoam  insulation.  On  the  downstream  side  of  the 
holder  is  a  plate  containing  the  thermocouples  used  to  measure  t^.  The 
test  section  into  which  the  holder  slides  contains  the  heaters,  pressure 
taps,  and  the  thermocouple  set  downstream  from  the  heaters,  which  measures 

V 

Inlet  Cone  and  Flow  Straightener 

This  section  was  designed  and  tested  by  Piersall  [15],  and  it  pro¬ 
vided  a  uniform  velocity  profile  to  the  air  entering  the  matrix. 

Pressure  Measuring  System 

Pressure  taps  are  located  in  the  test  section  upstream  and  down¬ 
stream  from  the  matrix  holder,  (see  Figure  5).  Two  other  taps  are  located 
at  the  orifice  section.  Each  pressure  tap  is  connected  by  flexible  tub¬ 
ing  to  its  corresponding  manometer  and/or  draft  gage.  The  following  in¬ 
struments  were  used  interchangeably,  depending  upon  the  flow  rate: 

1.  Ellison  Draft  Gage  Company,  0-3  in.  inclined  gage 

2.  Ellison  Draft  Gage  Company,  10  in.  manometer 

3.  Ellison  Draft  Gage  Company,  20  in.  manometer 

4.  Merriman  Instrument  Company,  120  in.  manometer 

5.  Precision  Thermometer  and  Instrument  Company,  mercury 
barometer 

Temperature  Measuring  System 

Temperatures  in  the  system  are  measured  at  four  locations:  (See 
Figure  5) inlet  to  the  system  (t^),  downstream  from  the  heaters  (t^). 
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downstream  from  the  matrix  it  ),  and  at  the  orifice  (t  ). 

3  0 

Temperature  t^  is  measured  by  two  different  sets  of  30  gage  iron- 
constantan  thermocouples,  which  were  made  by  Traister  [16],  Each  set 
consists  of  5  thermocouples  connected  in  series.  Each  thermocouple  is 
individually  wrapped  in  teflon  tape  to  prevent  shorting.  All  ten  are 
contained  in  an  open-faced  aluminum  tube  mounted  in  a  frame  at  the  exit 
of  the  inlet  cone.  The  aluminum  tube  shielded  against  radiation  from  the 
heaters. 


Temperature  t^  was  measured  by  a  set  of  five  .001  inch  diameter 
iron-cons tantan  thermocouples  connected  in  series.  The  output  of  these 
waa  bucked  against  one  of  the  sets  measuring  t^  so  that  the  output  of  the 
two  aeta  measured  (t2~t^). 

Temperature  t^  was  measured  in  the  same  manner  as  t^  and  the  output 
bucked  against  the  ether  set  measuring  t,  so  that  the  output  measured 

4> 

The  outputs  (t^-tj)  and  (t^-t^)  could  be  recorded  separately  or 
bucked  against  each  other  to  give  (t^-t^.  The  desired  outputs  were  then 
led  into  an  Aatrodata  Model  886  Wideband  Differential  D.  C.  Amplifier, 
where  they  were  amplified  100:1.  From  there  the  signal  was  led  into  a 
Brush  Mark  280  Strip  Chart  Recorder. 

Temperature  t^  waa  measured  by  a  30  gage  copper-constantan  thermo¬ 
couple  referenced  to  an  ice  junction.  The  output  was  read  on  a  Leeds  & 
Morthrup  Millivolt  Potentiometer. 

Heater  Power. 

The  power  to  the  heaters  was  supplied  from  the  250  V.  D.  C.  source 
in  the  laboratory.  The  actual  voltage  was  controlled  both  at  the  supply 


panel  and  at  the  apparatus  by  means  of  a  50  JL.  wound  resistor. 
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APPENDIX  B 


F  *  thermal  expansion  factor 
a 

)'  "  P/RT  the  specific  weight  of  fluid  flowing,  assuming  a 
perfect  gas  'lbf/f  t3) 

Y  ®  expansion  factor 

F  *  absolute  static  pressure  at  orifice  (Ibf/sq  ft) 

R  *  gas  constant  for  air:  53.35  (ft-lbf/lbra°R'* 

T  *  absolute  temper lture  at  orifice  (deg  R) 

Afp  “■  pressure  drop  across  the  orifice  in  inches  H^O 
Substituting  the  expressions  ter  K  and  in  equation  (B-5)  yields: 


m  =  35-s  _JL 

■/ t )-/3 

From  [1]  Fig.  40A 


&  £  /  I  ARP 


(B-6) 


y  -  /  —  (b.v/  -f-  a  35  )  _X. 

k  »  1,4  for  air,  ratio  jC-f>  j JRv- 

Also,  from  [1],  Fig.  38 

F  *  1.0 
a 

P  -  (P  -P  /13. 6) <0.4912*144)  lbf/ft2 
atm  o'* 

P  „  *  local  atm08*'heric  pressure  in  inches  Hg 

atm  r  ° 

Pq  *  static  pressure  upstream  of  the  orifice  plate  in  inches  Hg 
Making  the  above  substitution  in  (B-6)  with  he  constants  necessary  to  be 
dimensionally  consistent  yields: 

I— 

\  /-(?/ 


-i  ,035 


■sfiA 

•  A 


Aft 


-p  -  p 

(  txim  ° 


£) 

/  w 


Wot 


®/i3.d 


t„  f  45-9,7 


l/z 


(B-7) 
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Reynolds  Number 


Reynolds  Number  is  defined  as: 

Afe  =  A/  & 

A 

where  G,  the  mass  flow  velocity  *  m/Ac  *  m/pA^ 
and  Jj.  is  the  Liuid  viscosity. 

Substituting: 

Pc  =  rh  dh  /a  R+,  P  -  R-Ro/pi  fl,r  & 


(B-8) 


(B-9) 


M  =  tHnL./ Op. 

Maximum  Slope 


(B-10) 


The  maximum  slope  of  the  downstream  cooling  curve  (t^-t^)  iu  a  unique 
function  of  N  [13]  and  ^[6]. 


dfe=k)  =  0  .  * ) 

.  «<-(Va44  I  max 


^  ,  Ntu  a..d  ^  have  been  previously  defined,  but  for  convenience  are 


restated : 


q  }  M  =  ^  A 


Ws<i 


trur. 


rhtfL 


the  new  temperature  introduced  is  t^  and  is  equal  to  the  temperature 
of  both  the  solid  and  the  fluid  prior  to  cooling. 

Therefore : 


2"  /A(l^  ~  OvyCf  / Ws.Cs)  & 
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and 


cL  =  (Mb/  )ole 


letting: 


and 


*  mc^  (Fluid  stream  thermal  capacity  rate) 


C  *  W  c  (Matrix  heat  capacity) 

3  S3 


equation  (B-ll)  becomes „ 

oi(*/A^)  =  (CJCS) 


Furthermore 

^3  -  tc 


_  i:i  -  i  i 


*'  -<i  i,  -t; 

whose  derivative  is: 


4  / 


,  /  ^  6s  j 

’  £  3-t- 


combining  equations  (B-ll)  and  (B-12)  yields: 


cLCt/Ak.) 


MAX 


d( -t. ) 
oLo 


(B-ll) 


(B-12) 


(B-13) 
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From  the  above  sketch  of  the  downstream  temperature  response  curve; 

d-  <Vt)  , 

*  Y/X 

d  Q 


cLq 

\!  chart  speed 


d(t3-tx) 

ti~tl 

combining  with  the  ratio  oft 


matrix  capacity 


-  Y 

-  A  t 


C 


flow  stream  capacity  rate 
and  equation  (B-13) 


s  (l/8ec) 


4z3t) 

d(2/M) 


—  f  'Y  Char{  SpeecL 

£ f  At  X 


<B  14) 


This  value  of  maximum  slope  and  /\  are  then  used  to  enter  Table  III 
or  Figure  14  or  15  to  get  the  corresponding  value  of  N^. 

(2)  Centroidal  Technique 

For  the  centroidal  technique  the  recording  trace  of  (t^-t^ 
used.  Two  methods  are  actually  involved  in  utilizing  this  trace.  One 
method  i s  to  copy  the  trace  physically  with  carbon  paper  on  to  a  thin 
piece  of  cardboard.  Then  the  trace  copy  is  cut  out,  and  by  means  of 
a  "plumb-line"  attached  to  a  pivot,  which  allows  the  trace  to  swing  free, 
the  line  passing  through  the  pivot  point  and  the  center  of  mass  is  determined. 
This  is  done  for  several  points  (at  leant  3)  and  the  intersection  of  the 
different  lines  is  the  centroid  of  the  area  under  the  response  curve. 
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These  lines  are  shown  in  Figure  16.  Ideally  all  lines  will  intersect  at 
the  same  place,  however,  this  is  not  the  case  in  actuality  due  to  errors 
in  manipulation. 

An  alternate  method  is  also  used  to  find  'he  centroid  of  the  function 
(t^-tj).  This  involves  the  use  of  a  computer  program,  which  when  fed 
actual  curve  data  (normalized)  computes  the  centroid.  This  is  done  in 
the  absence  of  an.  apparatus  which  could  convert  directly  the  temperature 
response  to  information  useable  by  the  digital  computer,  such  as  a  paper 
tape  puncher. 

Once i*  determined  from  the  trace  of  <.t-j-t2)  it  is  necessary 
to  compute  1(g)  in  the  manner  described  in  Appendix  C.  This  value  of  1(g) 
is  then  used  to  determine  JX0  .  Recall  Equation  (14) 

/a™  =  A  +  .  99  1(f)  i  . O / 

or 

~  *99  Z(^)  +  .01 

Once  y-to  is  determined,  Kohlmayr's  curve  of  Nfcu  vB.yi^  ,  Figure  1, 
may  then  be  entered  to  give  Ntu» 

(3)  Zero  Intercept  Technique. 

To  interpret  the  physical  data,  i.e.,  the  recorded  trace  of 
(t^-tj),  which  represents  t^,  in  the  manner  described  in  Section  3,  of 
this  report: 

(i)  Extrapolate  by  means  of  a  French  curve  the  response 
curve  (t^-tj)  from  some  position  on  the  trace  where  the  transient  response 
of  the  heaters  has  died  out,  back  to  the  vertical  line  passing  through  the 
point  y*-  »  0.  This  point  would  be  the  position  where  the  first  change  in 
temperature  was  detected  by  the  thermocouples  and  is  referred  to  as  (t^-t^)0 
or  tf  (Ntu  ,0  >,  Figure  2B. 

(ii)  Record  the  value  of  (t^-t^)  at  the  intere^ tion  of 
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the  time  zero  line  and  the  extrapolated  response  curve,  also  record  the 


value  of  (t^-t^)  at  the  time  just  prior  to  the  heaters  being  deenergerized. 
This  temperature  is  referred  to  as  (t^-t^)^^. 

(iii)  N  is  then  equal  to, 


-ln(l-  <t3-fc1>0/ <t3-t1)inax> 


which  was  developed  in  Section  3,  as, 


-  Jbn.  'Mu,0)) 
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APPENDIX  C 


To  determine  the  transient  response  of  the  N?S  facility  for  differ¬ 
ent  heater  sets,  runs  were  conducted  in  which  the  temperature  difference 
(tj-t^)  was  recorded  as  a  function  of  time.  This  temperature  difference, 
which  was  previously  called  ’’upstream  temperature"  is  the  difference  in 
temperature  of  the  air  entering  the  test  rig  and  the  air  leaving  the 
heaters. 

Analytical  determination  of  the  heater  time  constant  is  straight¬ 
forward  and  presents  no  difficulty  based  strictly  on  the  analysis  of  a 
cylinder  in  cross  flow.  The  analytical  determination  of  the  thermocouple 
time  constant  is  more  difficult  and  requires  some  judgment.  Since  the 
thermocouples  are  constructed  of  .001  inch  diameter  wire,  it  is  not  pos¬ 
sible  to  describe  precisely  and  mathematically  the  geometry  of  their 
Junction.  In  the  construction  of  the  thermocouples  the  two  different 
wires  used  crossed,  at  an  angl«  approaching  180  degrees,  and  arc-welded 
at  the  junction.  Therefore,  the  junction  was  neither  spherical  nor 
cylindrical  (which  is  the  case  when  the  wires  are  butt-welded).  For  this 
reason  the  model  of  the  junction  was  arbitrarily  chosen  to  be  spherical 
and  calculations  were  based  on  a  spherical  junction  .002  inches  in  dia¬ 
meter.  On  this  basis  it  was  determined  that  the  time  constant  of  the  .003 
inch  heaters  was  nearly  three  times  greater  than  the  time  constant  of  the 
thermocouples,  while  the  time  constant  of  the  .001  inch  heaters  was  ap¬ 
proximately  two- thirds  that  of  the  thermocouples. 

Since  the  response  of  the  heater-thermocouple  circuit  wou^d  ^e  the 
sum  of  at  least  two  exponential  terms,  which  would  be  difficult  to 
separate  experimentally,  the  terms  have  been  treated  here  as  one  expon¬ 
ential  term,  which  is  a  good  approximation  of  the  actual  system  response 
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if  one  neglects  the  initial  points  on  the  response  curve.  From  the 
normalized  temperature  response  curves  the  following  empirical  relation* 
were  developed: 


tz  ~t,  _ 


where  Q  *  time  in  seconds  and 

is  the  heater  time  constant 
in  seconds  (C-l) 


/oj 

(t  1 


-  _  Q 


(C-2) 


where  the  value  of  — : - is  the  reciprocal  of 

log  r  (t2-t1) 

(t~*t^)raax. 

- 

the  slope  of  the  temperature  response  curve. 

The  heater  time  constant  was  then  plotted  against  mass  flow  rate  ( ) . 
From  these  curves  the  following  relationships  were  obtained: 

log  On  -  -  a  1  7,/y,  +  log  c  where  a  c  are  (C-4) 

G„  -  JCr>7~a'  constants  (C-5) 

Since  the  centroid  technique  makes  use  of  the  deviation  from  step  defined 
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Therefore  to  determine  1(g)  for  &  particular  run,  enter  the  curve  of 
vs.  /V)  (Figure  12)  with  apprcpri  .te  hh  .  Note  that  this  curve  is 
independent  of  the  matrix  used  in  the  run.  Then  multiply  by  the 

appropriate  ©<  to  obtain  1(g)  used  in  centroid  technique.  Figure  12 
shows  the  experimental  values  of  vs.  />?  plus  the  theoretical 

plots  of  heater  and  thermocouple  transient  response  for  both  sets  of 
heaters. 
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Figure  5.  Schematic  of  Test  Apparatus 
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Best  Available  Copy 


Core  T20-38  Test  Results 
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Figure  10,  Core  Solar  Ij.  Teat  Reaultn 
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Figure  15.  N^u  Versus  Maximum  Slope  With  Conduction 
Parameter  (Max.  Slope  .3-1.7) 


71 


Area  of  uncertainty 


Figure  18.  Photo  of  .003  In.  Wire  and  ,001  In.  Wire  Heater  Frames 
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